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We study the quantum mechanical correlation between two identical neutralinos in the decays of 
minimal supersymmetric standard model (MSSM) scalar tau (stau) pair produced in e"'"e“ anni¬ 
hilation. Generally, the decay products of scalar (spinless) particles are not correlated. We show 
that a correlation between two neutralinos appears near pair production threshold, due to a finite 
stau width and mixing of the staus and/or neutralinos, and because the neutralinos are Majorana. 
Because the correlation is significant only in a specific kinematical configuration, it can be observed 
only in supersymmetric models where the neutralino momenta can be kinematically reconstructed, 
such as in models with i?-parity violation. 
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I. INTRODUCTION 

Neutralinos are mixtures of the spin-1/2 superpartners 
of the neutral gauge and Higgs bosons in supersymmetric 
extensions to the Standard Model. The lightest one (x°) 
is generally considered to be the lightest supersymmetric 
particle (LSP), which is stable in i?-parity conserving sce¬ 
narios. Signals of supersymmetry ^USY) are expected 
to be observable at hadron Hiii or future e+e 
colliders Hi. 

In the minimal supersymmetric Standard Model 
(MSSM), neutralinos are Majorana fermions whose an¬ 
tiparticles are themselves. If SUSY (and the associ¬ 
ated neutralinos) is observed in nature, the determina¬ 
tion of their properties under CP conjugation, Majorana 
vs Dirac, will be a key issue: it is intimately related to 
the dimension of SUSY, e.g. N = 1 or N >2 SUSY 0. 

Many studies of the Majorana nature of neutralinos 
have been performed Reference H observed 

that the CP and CPT properties of the final state deter¬ 
mine the Majorana nature of neutralinos in the process 
^ XiXi ^ (* = 2,3,4). Selectron pair 

production in e“e“ scattering, e“e“ ^ e~e~ , would also 
be of great utility, as it could proceed only on account 
of t-channel Majorana neutralino exchange M- In ad¬ 
dition, Majorana effects in two- or three-body decays of 
neutralinos are known to exist M- 


‘Electronic address: kaoru.hagiwara@kek.jp 
^Electronic address: kentarou@post.kek.jp 
t Electronic address: rain@pas.rochester.edu 
5 Electronic address: tstelzer@uiuc.edu 


In this article, we study the Majorana nature of neu¬ 
tralinos in lighter scalar tau (stau) pair production and 
subsequent decay into r lepton (tau) plus LSP neutralino 
in e“'"e“ annihilation. 


e^e 


fi+fi 


r+T-x?X? , 
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via the quantum mechanical correlation between two 
identical neutralinos which exists only if they are Ma¬ 
jorana fermions. That a correlation should exist should 
be evident simply due to proper treatment of the Fermi 
statistics of the final-state neutralinos, but it is always ig¬ 
nored: conventional wisdom says that the correlation be¬ 
tween decay products of scalar (spinless) particles, such 
as staus, is absent. Because of this, little attention in the 
literature has been given to the expected Majorana effect, 
even though the process we consider here is among the 
most promising to study SUSY particles and many stud¬ 
ies of it have been performed [12, [l2 for exploration of 
other aspects. Naturally, no correlation is possible once 
a spinless particle is separated from the other particles 
at a macroscopic space-time distance. In other words, 
any interference effect disappears in the limit of long life¬ 
time of the spinless particles. We should hence consider 
the finite stau width, and investigate the region near pair 
production threshold. 

In addition to the finite width, the left-right mixing of 
the staus and/or the gaugino-higgsino mixing of the neu¬ 
tralinos are necessary for the correlation to be significant. 
This is because in the massless neutralino limit, same- 
helicity neutralinos are produced only when the taus have 
the same helicity, which occurs only when the fp-Tn mix¬ 
ing and/or the neutralino mixing are significant. We will 
discuss this in detail in the following section, and simply 
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note in passing that significant interference effects are 
expected in and production processes, 

which quantitative study will be reported elsewhere. 

The article is organized as follows. Section m gives a 
brief introduction of the mass eigenstates of the staus 
(fi, 2 ) and the gaugino-higgsino mixing of the neutrali- 
nos. The coupling is also given. In Sec. Iffll we 

present scattering amplitudes for the process e“'"e“ —> 
—> T'^T~XiXi and define the kinematical variables 
relevant to our analysis. Section CYI gives the total cross 
sections around the fi pair production threshold. In 
Sec. El we study in detail the kinematical correlations 
due to identical-particle interference effects and compare 
them to the Dirac (noninterference) case. Section Ivll is 
devoted to a summary and discussions. 

II. STAU AND NEUTRALINO MIXING 

In this section, we briefly introduce the staus’ left-right 
mixing and the neutralinos’ gaugino-higgsino mixing. We 
also give the ti-t-Xi coupling, where fi is defined to be 
the lighter stau mass eigenstate after L-R mixing. 

Because of the non-negligible tau Yukawa coupling, 
mixing occurs between the weak eigenstates tl and tr 
to form mass eigenstates fi and T 2 (rufi < rrif^) as 

(fL\ _ ( cosOf sinOA frA , . 

[jrJ y-sin6»f cos8fJ \T 2 J ' 


a± are expressed in terms of the f^-T-Ai and th-t-Ai 
couplings, and a^, respectively, as 

a± = cos 6 f ■ -I- sin0f • . (7) 

Here, 

a-=-^{U 2 i + Uiita,n 9w), a^ = U^iYr, 
a+ =-^2 C/fi tanOw, = U;^ Y, , ( 8 ) 

where Yr = — grriT/{V^mw cos P) and g is the SU{2)l 
gauge coupling. The couplings ab. and come from the 
stau-gaugino interactions, where chirality is conserved. 
On the other hand, and are proportional to the 
tau mass, and come from the f-r-H^ Yukawa coupling, 
which flips chirality. For further details see Ref. |lill| . 

At this point it is worth noting that either the left- 
right stau mixing or the gaugino-higgsino mixing of the 
neutralinos is needed for the Majorana interference ef¬ 
fects to be significant. This may be explained as follows: 
if both mixings are absent, the tau produced in decay is 
either left-handed (gaugino interaction) or right-handed 
(higgsino interaction), and the r’*' is of the opposite chi¬ 
rality. The two neutralinos then have opposite chiral¬ 
ity and their interference vanishes in the massless limit, 
where helicity coincides with chirality. 

III. AMPLITUDES AND KINEMATICS 


Neutralinos are mass eigenstates of the neutral gaugi- 
nos, B and H 3 , and the neutral higgsinos, 77° and 77° . 
The mass matrix in the X = {B, W 3 , 77°) basis is 

/ Ml 0 —mzswCfj mzSwSf3 \ 

^ ^ 0 M 2 mzcwC0 -mzcwS/B I 

-mzswcg mzcwcg 0 -g I 

V mzswsfj -mzcwsi3 -g 0 / 

(3) 

with the abbreviations sw = sin^W) cw = cos 9w, sp = 
sin/3, and cp = cos/3. Here, Mi and M 2 are the gaugino 
masses, g is the higgsino mass, and tan/3 = (77°)/(77°) 
is the ratio of the vacuum expectation values of the two 
Higgs doublets. The mass eigenstates are given by 

X° , (4) 

where U diagonalizes the above mass matrix as 

U'^M U = diag(TO^o, , m^o, m^o) . (5) 

Here, 0 < m^o < m^o < m^o < m^o . 

' — A 1 — A 2 — A .3 — A 4 

The decay ti —> tx° depends on both the left-right 
stau mixing (9f) and the neutralino mixing (Uij). The 
Ti-T-x° coupling can be expressed as 

£ = {a-P- + a+P+) pr + h.c., (6) 

with the chiral-projection operators P± = 4(1 ± 75 ). We 
denote left-handed (L) by ’ and right-handed (7?) by 
‘-f’ for notational convenience. The complex couplings 


In this section we present helicity amplitude formulae 
for the process 

e“(fc. A) -I- e'^(k, A) 

^ rr + A 

T“(fci,Ai) -f T+(fc2,A2) -f x?(pi,cri) +Xl(P2,a2), 

(9) 

where the four-momentum and helicity of each particle 
are defined in the center-of-mass (CM) frame of the e+e” 
collision. If neutralinos are Majorana fermions, the two 
neutralinos (x°) in the final state are identical. Therefore 
the crossed diagram (b) of Fig. E should be added to 
diagram (a) before the amplitude is squared. The relative 
sign between these two diagrams appears due to Fermi 
statistics. 



FIG. 1: Feynman diagrams for the process e''"e —> ^ 
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The full amplitude can be expressed as the product of the stau pair production amplitude {Mf), the two Breit- 
Wigner propagators for the staus, 

, ( 10 ) 

and two fi ^ decay amplitudes and That is, 

M =A 1 (A; Ai, A2, CTi, 172) 

= Mf{X;qi - q2)Df{qj )Df{ql ) M^{qi ; Xi,cri) M^{q2 ;\2,(T2) 

- Mf{X; q[ - 4 ) Df{q[^)Df{q2^) M^{q[-, Xi, a2) M^{q2; A2, tri) 

=Mi-M2 . (11) 


The intermediate stau momenta can be written in terms of the final-state particle momenta: qi = ki -fpi, (72 = ^2 +P 2 , 
q[ = ki + p 2 , and (72 = ^2 + Pi. The stau pair production amplitude is given by 


Alf (A; q) — — 1 H- 

s L s — m 


2 -^—-— PA {g- cos^ 9f + g+ sin^ 9f) v{k, - A) fiu{k, A) 

+ imzL z J 


( 12 ) 


with Z boson couplings to left- and right-handed charged leptons, g- = 9+ ~ respectively. Here, we 

neglect the electron mass and take A = —A. Using the straightforward Feynman rules for Majorana fermions given in 
Ref. 0 , the stau decay amplitudes for Aii are written as 


Afx(9i; Ai,cri) = u(fci, Ai) {a*_P+ + a+P_) w(pi, cri), 

Mx[q2-,X2,(J2 ) =u{p2 ,(T2) (a-P- + a+P+)v{k2,X2). (13) 

Similarly, for Ad 2 

M^{q[-, Ai,(T 2 ) = M(fci, Ai) {a*_P+ + a\P-)v(p 2 ,U 2 ), 

Ad X ( 172 ;''' 2 , 0-1 ) =u{pi,(Ji){a-P- + a+P+)v{k 2 ,X 2 ) ■ (14) 


In this article, we assume the lighter stau (fi) to be the next-to-lightest supersymmetric particle (NLSP), a common 
occurrence in many MSSM scenarios. Hence, all fi’s decay into r plus x?, and the total decay width Ff^ is just the 
partial width r(fi ^ tXi)- Using the fi decay amplitude Ad^^ in Eq. (ICT . the decay width is given by 

Tfi = r(n ^ Tx?) = y ^ |Adx(gi;Ai,cri)pd $2 = — (|a_p-f |a+p) mfi^l - , (15) 

Ai,cri '’’1 


with rUr = 0 and the Lorentz-invariant phase space factor 


(p = Yl n 

i—1 i—1 i—1 


d? ki 


(16) 


Throughout our study we neglect the tau mass, except in Yr of the fi-r-Xi coupling, given in Eq. 

Let us now define the kinematical variables. In the CM frame of the e“'"e“ annihilation, we choose the fi momentum 
direction as the z-axis, 

91 = ^(l + ^,0,0,/3), 

92 = #(1 + 2 ^, 0 , 0 ,-/ 3 ), ( 17 ) 

— 2 2 _ 

where /3 = /3(^, ^) with /3(a, 6) = (1 -|- -I- 6 ^ — 2a — 26 — 2 a 6 )^A^ and we choose the k x qi direction as the y-axis. 

For computational convenience, we parametrize the momenta of t~ and Xi with pi in the rest frame of 91 , 


k* _ sin 0 * cos(/)i,/3i sin 0 * sin(/))(,/3i cos^*) , 

Pi = (1 + ~/3r sin^C cos(/)i, —/3i sinS)” sim/i*, —/3i cos^*) , 


(18) 
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with = 1 — m-o/( 7 i. Those of and Xi with p 2 are then in the q 2 rest frame, 

^2 = ^ (1 - ^>/^2 sin6»* cos(/)^,/3Jsin6»* sin(/>*,/?* cos6»*) , 

P2 = ^ (1 + ^>-/32 sin6i2COS(/)2,-/32 sin6»2 sin(^2>-/32 cos6»2) , (19) 


with /JJ = 1 — miolqi- The two frames differ only by 

Xi 

a boost along the z-axis (see Fig. El- All the frame- 
dependent variables with a star superscript (*) are those 
in the fi rest frame. 



FIG. 2: Schematic view of the coordinate system. 

Before turning to the numerical study, it is worth¬ 
while to mention the relation of the helicities between the 
taus (Ai_ 2) and the neutralinos (cri,2). Generally, when a 
scalar (spinless) particle decays into two fermions, they 
always have the same helicity in the rest frame of the par¬ 
ent, due to helicity conservation in gauge interactions. In 
the massless neutralino limit, this relation remains even 
in the e+e” CM frame since the spin cannot flip due to 
the boost, i.e., 

Ai = CTi and A2 = cr2 for A 4 i, 

Ai = (72 and A2 = CTi for A42 ■ (20) 

Therefore, in the massless neutralino and tau limit, when 
T~ and T+ have the same helicity (Ai = A2), the two neu¬ 
tralinos also have the same helicity, and we expect sig¬ 
nificant interference effects between the two amplitudes, 
A 4 i and A42. On the other hand, there is no interference 
for the Ai ^ X2 case since either the amplitude A 4 i or 
A42 vanishes. For finite neutralino mass, however, the 
simple relations between the helicities in Eq. do not 
hold due to the spin-flip effects by Lorentz boosts. In¬ 
terference effects can then appear even for the Ai ^ X2 
case. In the following, we always sum over neutralino 
helicities, but not tau helicities. Notice that the tau po¬ 
larizations are observable statistically from their decay 
distributions. 

In order to confirm the above kinematical analysis, we 
consider the interference term analytically. In the spin- 
summed squared amplitude = Sl-^i ~-Ad2p 

(cf. Eq. lUTll b the interference term is given by the real 
part of M1M2 as 

/ = -2Re M1M2, (21) 

A,(Ti ,(72 


where we sum over the initial electron polarizations (A) 
and the neutralino polarizations (cti and CT2), but keep 
the helicities (Ai and A2) fixed. The minus sign is the 
result of Fermi statistics. The interference term I can be 
expressed as the product of the production part Ip, the 
Breit-Wigner propagator part, and the decay part Ip- 

/ = 2 Re (Ip • D,{ql)D,{ql)D%{q[^)D%{q'Y ■ Ip) , 

( 22 ) 


where 


Ip=-Y ^riX; qi - <Z2) M;(A; q[ - g'), ( 23 ) 

A 

Ip = ■Adx((?i;Ai,cTi) Alp^(g 2 ;A2,(72) 

X AIJ(g(; Ai,( 72 ) A 4 ^(g 2 ; A 2 ,( 7 i) . ( 24 ) 

Note that we include the minus sign explicitly in the pro¬ 
duction part. For the decay part, using the fi decay 
amplitudes of Eqs. and da, the relations 

^■u^(p,s)'u(p, s) =( 7 "^(^-to), 

S 

^u(p, s)m'^(p, s) = (^-to)( 7 '^, ( 25 ) 

S 

and the properties of the charge conjugation matrix C, 
we obtain 

_ J|a-P|a+P tr[/^i^i )(f2^2] for Ai = A2 , 

^ [5 w|o (|a-|^ + |a+|'‘) tr[ |iii/^2] for Ai 7^ A2 . 

( 26 ) 


This confirms the above kinematical analysis. 


IV. TOTAL CROSS SECTIONS 

In this section, we present total cross section for the 
process e“'"e“ —*■ around and well below the 

fi pair production threshold. In addition to the double¬ 
resonance contribution (Fig. ^ which we have discussed 
so far, we further consider the single-resonance contribu¬ 
tions to the final states, shown in Fig. El There are two 
additional diagrams, each having its own crossed diagram 
for the Majorana neutralino case. 

The cross section for the process e+e” ^ t~^t~XiXi 
averaged over initial electron polarization and summed 
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FIG. 3: The single-resonance contributions to the final state 
in e'^'e” annihilation. 


over neutralino polarizations reads 



X.<7i ,(72 

(27) 

where Ai and A 2 are the helicities of r“ and r+. In 
addition to the initial-state spin-average factor, we di¬ 
vide by the statistical factor for two identical neutralinos 
in the hnal state. The amplitude A4 is summed over 
all diagrams; not only the double-resonance contribu¬ 
tions in Eq. d, but also the single-resonance contri¬ 
butions. The four-body phase space factor d ^4 can be 
decomposed as the two-body phase space d^ 2 , and can 
be parametrized by the kinematical variables defined in 
Eqs. (|T7 |l - P)l as 


d$4 =c?<i>2(fc -l- k 

X d^2{qi 

27r(327r^)3 


= 9i 


dqf dql 


= ki +pi)d(^2{q2 = k2 +P2) 

dql dql d cos 9 d cos 9\dtjil d cos 9^ d(j )2 , 

(28) 


where 9 is the scattering angle between the e momen¬ 
tum {k) and qi in the e+e“ CM frame (see Fig.|21l. 

Now let us explain several parameters we use in our 
numerical analysis. All the helicity amplitudes, includ¬ 
ing the single-resonance contributions, are calculated by 
HELAS subroutines |l fT|| . and numerical integrations are 
done with the help of the Monte Carlo integration pack¬ 
age BASES [ 3 . We fix the fi mass at rrif^ = 150 GeV. 
As for the SUSY parameters, including the left-right stau 
mixing, we take the following values so that the Majorana 
effects are expected to be large: 

9f = 45° , tan/3 = 50, M 2 = 300 GeV, p = 70 GeV, 

(29) 

and adopting the mass relation Mi = | tan^ 9 wM2. This 
parameter set corresponds to a higgsino-like neutralino 
LSP. For these parameters, the xl mass is = 50 GeV, 
which yields a n total decay width of Tf, = 0.52 GeV, 
using Eq. (O. 

Figure 01 shows the helicity-dependent total cross 
section, Eq. 03, for the process e+e ^ r+r xlxl as a 


FIG. 4: Total cross section as a function of collision energy, 
of the process e"*'e“ ^ T'^T~XiXiy for (^) Ai = A 2 and (b) 
Ai 7 ^ A 2 , where Ai ,2 is the helicity of the Solid and dashed 
lines show the case that neutralinos are Majorana and Dirac 
fermions, respectively. Dotted lines are for each contribution 
from the ri pair and singly resonant fi production without 
the interference term. 


function of the GM energy (\/s) in e+e” collisions. In or¬ 
der to show the single-resonance contributions, the cross 
section is shown starting from rather low values. Solid 
lines denote the cross sections including the crossed di¬ 
agrams, such as Fig. Mb), which should be present for 
Majorana neutralinos. The dashed lines are obtained 
by neglecting the crossed diagrams, which corresponds 
to Dirac neutralinos. As a reference, each contribution 
from the fi pair and the single-fi production (without 
the interference term) is shown by dotted lines. 

Above the fi pair production threshold (y^ > 2mf^ = 
300 GeV) the contribution from fi pair production 
(Fig. lU is dominant. Below pair threshold the single¬ 
resonance processes (Fig.^J contribute dominantly, even 
though the cross section becomes very small. The reduc¬ 
tion of the cross section by the interference effects for 
Majorana neutralinos can be seen below the ti pair and 
the single ti production thresholds. Unfortunately, the 
total cross section in the region below threshold where 
Majorana interference effects become important is not at 
a level which could be observed. 


V. MAJORANA EFFECTS 

We have seen that it is difficult to observe a Majo¬ 
rana interference effect in the total cross section. In this 
section, we therefore study in detail the kinematical cor¬ 
relations due to interference effects that appear only for 
Majorana neutralinos. We present several distributions 
for the process e+e“ ^ t~^t~XiXi near the stau pair 
production threshold, and discuss the Majorana effects 
as a function of the finite ti width. We consider the 
following three distributions to see the interference: 

(a) cos 9*: defined in the qi rest frame in Eq. (Hl^ (See 
Fig. 13; 
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(fb) (fb/GeV) (fb/GeV) 



FIG. 5: Distributions of (a) cos^i, (b) Er, and (c) M{xx) for the process e'^'e” ^ XiXi fo \/s = 302 GeV 

for 1, 10,20 X Ffj in the massless neutralino limit. Solid and dashed lines show the Majorana and Dirac neutralino cases, 
respectively. Also shown is the behavior of the interference term by dotted lines. 


(b) Er'. the energy in the laboratory frame; 

(c) M{xx)- the invariant mass of the neutralino pair, 

= {Pi + P2Y = [k + k - ki - ^ 2 )^ . 

Note that, from the experimental point of view, these 
variables are not observables, since tans always decay 
into at least one neutrino, which escapes detection to¬ 
gether with the neutralinos; the kinematical system is un¬ 
constrained by observables and cannot be reconstructed. 
Therefore, our studies may be regarded as pedagogical, 
or may apply to models where the neutralino momenta 
can be kinematically reconstructed, such as those with 
i?-parity violation. 

To begin with, for simplicity we consider only the t\ 
pair production process e+e" —> 
and take the massless neutralino limit. 

Figure 0 shows the distributions of (a) cosd^, (b) Er, 
and (c) M{xx), at yG = 302 GeV for the Ai = A 2 case. 
We use the same SUSY parameter set in Eq. (123) as 
the total cross sections. To examine the finite fi width 
effect, we vary the total width Tfj^ as 1,10,20 x Tf^. 
In the limit of m^o = 0, the fi total decay width is 
Tfj = 0.66 GeV from Eq. (^3, hence 10 x Tf^ = 6.6 GeV 
and 20 X Tfj^ = 13 GeV. Above threshold, the wider 
the decay width, the larger the cross section. Solid and 
dashed lines show the Majorana and Dirac neutralino 
cases, respectively. Also shown as a reference is the be¬ 
havior of the interference term only, Eq. by dotted 
lines. For the realistic 1 x Ff^ case, the interference pat¬ 
terns are barely discernible in Fig. [3 The contribution 
of the interference term to the cross section is at most 
about 0.1%, 3% and 6% for 1,10 and 20 x Ff^, respec¬ 
tively. The interference effect is roughly proportional to 


the fi decay width. The following features are worth 
noting: (i) For the Ai yf A 2 case, there is no interference, 
as expected from the discussion of Sec. Hill (ii) For the 
Dirac neutralino case, the cross section does not depend 
on cos 91 or cos 9^ since the two produced staus decay in¬ 
dependently. On the other hand, for the Majorana case, 
the distribution is no longer flat. The cos 9^ distribution 
is the same as that of cos 9^, but for a relative sign, ac¬ 
cording to GP-invariance. (iii) A major portion of the 
interference pattern disappears when we integrate out 
the kinematical variables, such as cosSj", Er, or M{xx). 
This is why one can hardly view the correlation effect in 
the total cross section above pair production threshold 
in Fig. 21 (iv) As for the ^/s dependence of the interfer¬ 
ence effect, it exists even at higher energies. However, 
the cross section also grows, and the relative effect of the 
interference term becomes smaller. Therefore the Majo¬ 
rana effect can be seen only near threshold, even though 
the cross section is small due to the p-wave threshold 
factor of . 


Now we attempt to explain why the interference term 
changes sign as in Fig. |3 We take two different ap¬ 
proaches. One is to consider the interference term, 
Eq. (121 |l . analytically. Another is, more physically, to 
investigate them kinematically. 


First, let us consider the fi pair production part of the 
interference term, Eq. 123, analytically. Using the fi 
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pair production amplitude of Eq. HU), 


/p = - ^ ^ v{k, -A) (^1- ^2) A) 

^ xu{k,X){fi[-fi2)v(k,-X) 

= (30) 

The trace part is given by the dot-products of the kine- 
matical variables in Eqs. (|T7|l-ni?ll. The dependence on 
the scattering angle 9 appears only in this part. We in¬ 
tegrate it out and obtain 


IpdcosO 


4/36"^ 




cos 01 
cos 91 


- 


(31) 


where the neutralino mass is kept explicitly for later dis¬ 
cussion. This shows that the interference term behaves 
like — COS0*, or -|-cos 02 in the massless neutralino limit. 
This is consistent with Fig. [3 (a), where one can see neg¬ 
ative interference for cos 0 )“ > 0 and positive interference 
for cos 01 < 0 . 

Next, we turn to the kinematical approach to under¬ 
stand the Majorana effects more physically. The ques¬ 
tion is if there is a case when both amplitudes M.i 
and AI 2 of Eq. m are large, thus the total amplitude 
M = Ml—M 2 iS, suppressed due to the relative sign from 
Fermi statistics. For instance, the amplitude M should 
be suppressed when the two neutralinos have identical 
spins and four-momenta. The same spin is realized by 
setting the t~ and helicities equal in the massless 
neutralino limit, as previously noted. Let us determine 
the condition for which two neutralinos have the same 
four-momentum, i.e., pi = P 2 - It is obvious that the rel¬ 
ative azimuthal angle is zero, (j>* = — 4>i =0°. For 

simplicity, we consider the limit that both staus are on 
mass shell, qf = = Wfj, where the amplitude Mi is 

significant. In this limit, we can find the following simple 
kinematical point: 


cos 02 = —cos 01 , COS0*=/3, (32) 


with /3 = (1 — 4m|ys)^/^. At this point, not only q\ and 
but also q[^ and q'^ are on ti mass shell: 


2 2 / 2 / 2 2 

01 = 02 = 01 =02 = 


(33) 


We note that the propagator momenta squared of the 
crossed diagram are expressed as 


/ 2 
01,2 


(1 -|- /3^)(1 -I- cos 01 cos 02 ) 

-I- (1 — 0^) sin 01 sin 02 coscp* 


± 2/3(cos 0* -I- cos 02 ) 


(34) 



FIG. 6: Contour plot of the differential cross section, Eq. I45II . 
for the process e"''e“ —> r+fT —> T^r“xiXi a-t \/s = 302 GeV 
in the m^o = 0 limit when qf ~ q 2 = and (p* = 0°, is 
shown in the cos 0 i-cos 02 plane, where cos 0 i ,2 are the kine- 
matical variables defined in the gi ,2 rest frame. qi = rrif^ 
and q0 = are also shown by dashed lines. 


Figure El is a contour plot of the differential cross sec¬ 
tion for the process e+e“ ^ ^ ''"''’''■"XiXi at 

0s = 302 GeV in the massless neutralino limit when 
ql = ql = and (p* = 0 , 




dqidq 2 d COS 0i d cos 0^ d(p* 


(35) 


ql=ql=ml , 4>’-=q 


for the Ai = A 2 case. Using Eq. (El, we also show the 
00 = and q0 = trajectories by dashed lines. 
The interference effect can be seen along these lines, es¬ 
pecially around the intersection points, where both q0 
and q0 approach and the effect is largest. This arises 
from the double Breit-Wigner factor D{q0)D[q0) of the 
crossed amplitude M 2 - Note also that the sign of the in¬ 
terference term changes over the q0 = rri-^ 01 q0 = 
trajectory because of the Breit-Wigner resonant factor. 
Around one of the intersections (black region), which 
corresponds to the kinematical point of Eq. El with 
/3 = 0.115, the cross section is strongly suppressed be¬ 
cause the two neutralinos have the same momenta. On 
the other hand, it is enhanced around another intersec¬ 
tion (white region), i.e., (cos0)“, cos 02 ) = (—/?,/?). That 
is because the four-momenta of the t~ and leptons 
become identical at this point, and the two interfering 
amplitudes are constructive. Here one can also see nega¬ 
tive interference for cos 0 * > 0 (cos 0 ^ < 0) and construc¬ 
tive interference for cos 0 j‘ < 0 (cos 02 > 0 ). 

Let us now try to explain the consistency among the 
three distributions cos 0 j‘, Er and M{xx) in FigH In 
the phase space limit of M{xx) = 2 toxo = 0 , the 
four-momenta of the two neutralinos are identical {pi = 
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FIG. 7: Distributions of (a) cos0i, (b) Er, and (c) M{xx) for the process e''"e“ —> t'^t~XiXi for the Ai = A 2 case at 
= 302 GeV for 1,10,20 x Ff^. Solid and dashed lines show the Majorana and Dirac neutralino cases, respectively. Also 
shown is the behavior of the interference term by dotted lines. 


P 2 )- As discussed above, therefore, negative interference 
should be expected near this limit due to Fermi statis¬ 
tics. Furthermore, when = ( 72 ^ = rrif^ in this limit, 
it is obvious that the amplitude is suppressed in the re¬ 
gion COS 01 > 0 and cos 02 < 0. In addition, due to boost 
effects, this kinematical region corresponds to large Et- 

So far all our distributions have been given for massless 
neutralinos. Let us now show results for finite neutralino 
masses. 

Figure [71 is the same as Fig. O except we include the 
single-resonance contributions to the final state and a 
finite neutralino mass, m^o = 50 GeV. Since the single¬ 
resonance contributions are not so small just above the 
stau pair threshold (see Fig. O, the distributions are af¬ 
fected significantly. As for the M{xx) distribution, the 
lowest value of the invariant mass is 2m^o = 100 GeV in 

Al 

this case. We find that the interference pattern is basi¬ 
cally the same as in Fig. O where we consider only the 
Ti pair production process in the = 0 limit. 

However, one can also see a difference in the behavior 
of the interference between Fig.Qand Fig.O The region 
of constructive interference becomes larger than that in 
Fig. El The reason is given by Eq. eU: the third term 
proportional to the neutralino mass squared is additive 
with the first term, — cosSj", when we consider a finite 
neutralino mass. 

We can also repeat the kinematical analysis for the 
Majorana case, taking into account a finite neutralino 
mass. Equation is slightly modified by the boost 
effect as 

1 -I- mlojml 

cos 02 = — cos 9l , cos 9l = - 2 ^-;— 5 - P > P, (36) 

1 - 


where the two neutralinos have the same four- 
momentum. On the other hand, the condition that the 
momenta of t~ and are the same does not change, 
namely (cos0^,cos02) = {—PjP)- This supports the ten¬ 
dency of the interference pattern to increase with finite 
neutralino mass in Fig.|7| 

Figure|S|is the same as Fig.[7| but for the Ai A 2 case. 
As we discussed and explicitly showed in Eq. (ESI at the 
end of Sec. imi one can see the Majorana interference 
effect even for the Ai A 2 case. The distributions for 
the Ai = A 2 case are much more useful than those of the 
Ai y^ A 2 case because in (a) there is a significant shape 
change which includes two inflection points, whereas the 
distribution is monotonic in the latter case; and because 
in (b,c) the distributions for Ai = A 2 involve a peak shift, 
not just a magnitude change as for Ai y^ A 2 . 

A final question we may ask is, what would be the 
observability of the Majorana interference effect? We 
answer this by calculating the integrated luminosity re¬ 
quired at a future linear collider to observe a 3cr ef¬ 
fect, using the case MSSM described above and colli¬ 
sions at -^/i = 302 GeV. The simplest observable is 
the cos 01 distribution of Fig. 13a), which has a small 
forward-backward asymmetry A. We furthermore make 
the estimate using the nonphysical enhanced-effect case 
of 10 X Ffj. The formula for the statistical uncertainty 
on an asymmetry measurement is, 

AH = ^Nf{2Nb/N^)^ + Vb(2Vf/V2)2 , (37) 

where Nf{Nb) is the number of events with cos0* > 0 
(cos 01 < 0) and N = Np + Nb is the total number 
of events. Since the asymmetry is small, 0.02 in this 
case, we may use the approximation Np = Np = N/ 2 . 
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FIG. 8 : The same as Fig.Q but for the Ai 7 ^ A 2 case. 


Substituting and rearranging, we arrive at a formula for 
the required luminosity to observe a Scr effect: 

-^min — (3M)VfTtotal . (38) 

Plugging in = 0.02 and Ctotai = 1-2 fb, we arrive at an 
estimate of 18,750 fb“^. One is swift to conclude that 
even in an enhanced-effect scenario due to abnormally 
large finite stau width, a future linear collider unfortu¬ 
nately cannot observe this effect. 


VI. SUMMARY AND DISCUSSIONS 

We studied the quantum mechanical correlation be¬ 
tween two identical neutralinos, which exists only if they 
are Majorana particles, for the process e+e“ ^ 
T'^T~XiXi- We also considered the single-resonance con¬ 
tributions to the final state. 

We found that the correlation between two neutralinos 
appears near fi pair production threshold in the presence 
of a finite stau width and mixing of the staus and/or neu¬ 
tralinos. We discussed the finite width effect in detail, 
and found that the correlation effect tends to be propor¬ 
tional to the decay width. Distributions in several kine- 
matical variables, such as cos0*, Er, and M{xx) as de¬ 
fined in Sec. m show the interference effect, although the 
effect largely disappears after integrating over these dis¬ 
tributions. Because the correlation effects are significant 
only in a specific kinematical configuration, they can be 
observed only in models where the neutralino momenta 
can be kinematically reconstructed, such as in models 
with i?-parity violation. Unfortunately the interference 
pattern does not persist in the angular distribution of 
the final-state taus in the lab frame, disallowing these 


observables to be determined in an i?-parity-conserving 
MSSM scenario. 

Our brief estimate of the potential observability of 
the Majorana interference effect (assuming an i?-parity- 
violating scenario) using the asymmetry of the stau de¬ 
cay angular distribution in its rest frame is unfortunately 
not optimistic. It appears that an unrealistic amount 
of integrated luminosity at a future e“'"e“ collider oper¬ 
ating at stau pair threshold would be required. Thus 
this particular interference effect for stau NLSP pairs re¬ 
mains a pedagogical observation, but a very interesting 
one nonetheless. 

Before closing our discussions, we point out that a sig¬ 
nificant interference effect is also expected in and 

AiMfl PS'ir production processes. The quantitative study 
will be reported elsewhere. 
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